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Abstract: Using the generalized Konishi anomaly (GKA) equations, we derive the effective 
superpotential of four-dimensional M = 1 supersymmetric SU(A''c) gauge theory with Nf = 
Nc + 2 fundamental flavors. We find, however, that the GKA equations are only integrable 
in the Seiberg dual description of the theory, but not in the direct description of the theory. 
The failure of integrability in the direct, strongly coupled, description suggests the existence 
of non-perturbative corrections to the GKA equations. 



1. Introduction 



Holomorphicity of the supcrpotcntial and gauge couplings, global symmetries and the weak- 
coupling limit enable one to obtain exact results in supersymmetric gauge theories (for reviews 
see [1, 2]), making these theories more tractable than their non-supersymnictric cousins. 
Since supersymmetric gauge theories exhibit a wealth of non-perturbative phenomena such 
as dynamically generated superpotentials with associated confinement or chiral symmetry 
breaking [3, 4], deformed classical moduli spaces [5], Seiberg duality [6], etc., and since 
some of these phenomena also occur in non-supersymmetric theories, supersymmetric gauge 
theories are usually considered as a way to qualitatively study non-perturbative aspects of 
ordinary gauge theories. Therefore having a clear picture of the behavior of supersymmetric 
gauge theories may shed light on a better understanding of the dynamics of strongly-coupled 
gauge theories with no supersymmetry. 

Despite much progress in the effective dynamics of four-dimensional M = 1 supersym- 
metric QCD, the behavior of the effective superpotential for a number of flavors Nf large 
compared to the number of colors N^, is not well-understood. This is because, firstly, when 
the number of flavors increases there are typically additional light degrees of freedom at the 
origin of the moduli space that one needs to include in the effective description. Secondly, 
the effective superpotentials become singular when expressed in terms of the local gauge- 
invariant light degrees of freedom away from the origin; more precisely, the potentials derived 
from such effective superpotentials have cusp-like singularities at their minima [7]. Thirdly, 
the dependence of these effective superpotentials on the strong coupling scale of the theory A 
is such that they apparently diverge in the the weak coupling limit A ^ 0. Because of these 
problems the physical meaning of such superpotentials is thought to be problematic. 

We have argued elsewhere [7] that effective superpotentials for the light gauge-invariant 
degrees of freedom away from the origin must nevertheless exist. Furthermore, direct compu- 
tation in the case of SU(2) superQCD shows [7] that these superpotentials, although singular, 
are nevertheless physically sensible, and reproduce both the low energy physics as well as cer- 
tain higher-derivative terms in an intrinsic description on the moduli space away from the 
origin [8]. 

In this paper we extend our arguments to SU(iVc) superQCD by computing its singular 
effective superpotential. Unlike the SU(2) case, the SU(Arc) case has a smaller global sym- 
metry group, making it harder to find the superpotential. We deal with this by solving a 
system of differential equations for the effective superpotential [9] derived from the general- 
ized Konishi anomaly (GKA) equations [10]. The complexity of this system increases with 
the number of massless fundamental flavors Nf, but we are able to solve them in the first 
interesting case, Nf = Nc + 2. 

However, there are some subtleties involved in applying the GKA equations in this case: 
the GKA equations are not integrable when applied to the (strongly coupled) direct descrip- 
tion of the theory, but arc integrable when applied to the (weakly coupled) Seiberg dual 
description of the theory. In the rest of this section, we will explain these subtleties in more 
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detail and discuss the issues that they raise concerning the possible non-perturbative exact- 
ness of the GKA equations. We leave the technical details of the calculations to subsequent 
sections. 

The indirect argument for the existence of the effective superpotential referred to above 
goes as follows: Wilsonian effective supcrpotcntials are assured to exist only if there is a 
region in the configuration space of the chosen chiral fields where all of them are light to- 
gether and comprise all the light degrees of freedom. If this condition is satisfied, then the 
resulting effective superpotential can be extended over the whole configuration space by ana- 
lytic continuation using the holomorphicity of the superpotential. For a large enough number 
of flavors the theory becomes IR free and the only region where all the components of the 
chosen chiral vcvs become light at the same time is at the origin. We know what the light 
degrees of freedom are near that point since we have a weakly coupled lagrangian description 
there. The physics can be made arbitrarily weakly coupled simply by taking all scalar fleld 
vevs (0) << A where A is the strong coupling scale (or UV cutoff) of the IR free theory. In 
this limit the physics is just the classical Higgs mechanism, and all particles get masses of 
order {(p) or less. The Wilsonian effective description results from integrating out modes with 
energies greater than a cutoff, which we take to be some multiple of {(j)). The effective action 
will then include all local gauge-invariant operators made from the fundamental fields in the 
lagrangian and which can create particle states with masses below the cutoff. For the purpose 
of constructing the effective superpotential, the relevant local gauge-invariant operators are 
those in the chiral ring. It is then just a matter of constructing in the classical gauge theory a 
set of operators which generate the chiral ring. We will refer to this set as the classical chiral 
operators of the theory. 

In a weakly coupled SU(A'c) superQCD a basis of local gauge-invariant operators in the 
chiral ring (the classical chiral operators) is comprised of just the glueball, meson, and baryon 
operators [10, 11]. An effective superpotential which is a function of these operators must 
then exist. For Nf > Nc + 1, the quantum moduli space is also the same as the classical 
one [5], but effective superpotentials (singular or not) for these cases have not been found 
before. Also, this theory has an equivalent description in the IR in terms of a "Scibcrg dual" 
S\J{Nf — Nc) supersymmetric QCD with Nf (dual) fundamental quarks and anti-quarks and 
a set of singlet scalars coupled to the dual mesons through a superpotential [6]. 

Note that the above argument does not directly show the existence of such effective 
superpotentials in the asymptotically free case. In particular, for theories in the "conformal 
window" where neither the direct nor Seiberg dual description is IR free (^Nc < Nj < 3Nc 
for SU(A^c) gauge group), we have no useful description of the light degrees of freedom at 
the origin of moduli space. Nevertheless, given an effective superpotential for an IR free 
theory, one can then successively add mass terms to the effective superpotential and integrate 
out massive flavors to derive consistent effective superpotentials in the conformal window. 
This then assures us that effective superpotentials exist for all numbers of light flavors in 
supersymmetric QCD. 

Our method for deriving the effective superpotential for the Nf > Nc + 1 theory will be 
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to integrate the generalized Konishi anomaly (GKA) equations [10] following the approach of 
[12, 9]. The resulting equations become very complicated [9] for large numbers of flavors, so 
we are not able to solve them directly in the IR free case Nf > 3Nc, and then integrate out 
flavors as in the above argument. 

For Nf = Nc + 2, however, the GKA equations simplify to a first order matrix differential 
equation simple enough that we can analyze it. We show that the GKA equations for the 
effective superpotential are not integrable in this case for Nc > 2. This is not in direct conflict 
with the general arguments advanced above: for Nf = Nc + 2 and Nc > 2, the theory at 
the origin is strongly coupled in terms of its microscopic fields, so an effective description in 
terms of the chiral ring operators made from these fields simply need not exist. 

However, this failure of integrability presents a sharper puzzle in light of the following: we 
can nevertheless calculate an effective superpotential by using the fact that for Nf = Nc + 2 
with Nc> ^ the Seiberg dual description is an IR free SU(2) gauge theory. By applying the 
GKA equations to the Seiberg dual description we derive the effective superpotential of the 
theory in terms of the dual chiral fields. It is given in equation (4.18) below, where we have 
used the map [6] between direct and dual chiral operators to interpret this as an effective 
superpotential in terms of the classical chiral operators of the direct theory — the mesons Mj, 
baryons Bij and S*-', and the glueball S. 

We have thus found the effective superpotential in terms of the classical chiral operators. 
This raises the question of why were the GKA equations not integrable in the direct theory in 
the first place? We interpret this failure of integrability as indicating that the GKA equations 
get non-perturbative corrections. It remains to characterize more precisely the nature of the 
non-perturbative corrections to the GKA equations. One possibility is that there exists a non- 
perturbatively modified set of GKA equations in terms of the classical chiral operators (i.e., 
glueball, mesons, and baryons in our case). Another possibility is that there are additional 
operators in the chiral ring which are independent of the classical chiral operators; when 
included, they could render the GKA equations integrable, in much the same way that the 
extra singlet field in the Seiberg dual description does. It has not been ruled out that such 
additional fields could be seen in the semi-classical description as higher derivative chiral 
operators.^ (Note that the non-trivial higher-derivative chiral operators constructed in [8] are 
not candidates, since they are Q-exact when extended off the moduli space to the configuration 
space of chiral vevs.) Finally, both possibilities — a non-perturbative deformation of the GKA 
equations and the inclusion of additional chiral fields — could occur together. 

Once the issue of non-perturbative corrections to the GKA equations is raised, it applies 
equally well to the GKA equations derived in the Seiberg dual description. It is an open 
question whether the effective superpotential derived below in the dual description is correct 
or not for Nf > 6. For even though, when Nf > 6, the dual description is weakly coupled 
at the origin of moduli space, it becomes strongly coupled an arbitrarily small distance away 
from the origin since the superpotential term in the dual theory destabilizes the free fixed 

^We thank S. Hellerman for discussions on this point. 



-3- 



point at the origin [6]. 

The remainder of the paper carries out the computations described quahtatively above, 
and is organized as follows. To illustrate the method in a simple case first, and for later 
comparison to the Seiberg dual description of the SU(A''c) case, in section 2 we consider the 
case of SU(2) gauge group with Nf > 4 flavors. We show how the GKA equations written 
in terms of SU(2) mesons and baryons gives an effective superpotential matching that found 
in [7] where we worked instead with the single antisymmetric meson field appropriate to an 
Sp(l) description of the theory {i.e., one which makes the enlarged global symmetry group of 
the SU(2) compared to the general S'U{Nc) theory manifest). In section 3 we apply the GKA 
equations to SU(A^c) superQCD with Nf = Nc + 2 and Nc > 4. Wc show that the resulting 
equations for the effective superpotential arc not intcgrable for Nf > 6, but that they are 
integrable and match the SU(2) result of the previous section for Nf = 4. 

In section 4 we apply the GKA equations to the Seiberg dual of the theory in section 
3, and solve for the effective superpotential. The form of this superpotential is complicated: 
integrating out the heavy glueball gives an effective superpotential of the form y/detMf{X) 
where X = MBM'^B/detM, but a closed-form expression for f{X) is not found. Instead, 
we show that / obeys a nonlinear first order matrix differential equation (4.34). A power 
series expansion of the solution to order is computed in (4.29). We then compare this 
result to the SU(2) effective superpotential of section 2 when Nf = 4, and show that they 
agree at least to order X^. 

2. Weff for Nf > 4 SU(2) superQCD 

We show how to calculate the effective superpotential of SU(2) superQCD with Nf > 4 using 
the generalized Konishi anomaly equations, following [12, 9]. 

SU(2) superQCD has Nf massless quark and anti-quark chiral multiplets, Qa and Q^, 
transforming in the 2 and 2 of the gauge group, respectively. Here a = 1,2 is the color 
index and i = 1, . . . , Nf the flavor index. The apparent global symmetry of the theory is 
SU(A^/) X SU(A'^j) X U(l)s X U(1)r. Since 2 and 2 are equivalent representations, though, 
SU(2) superQCD actually has the larger symmetry group SU(2A^j) x which is large 

enough to determine the effective superpotential uniquely [7]. Here, since we are looking to 
the generalization to SU(A'c) which only has the smaller symmetry group, we will analyze the 
SU(2) case in terms of the Q's and Q's keeping only the SU(A^/) x SU(A^/) x U(1)b x V{1)r 
symmetry manifest. 

The classical moduli space is parameterized by the vevs of the meson M and the baryon 
B, B chiral superfields defined by 

M\ := Q\Q-^, B'^ := e^'QiQl, Bij := eabQfQ'j, S := tr(W^"T^«)/(327r2), (2.1) 

where we have also defined the glueball chiral superficld S. These fields can be assigned the 
charges R{S) = 2 and R{M) = R{B) = R{B) = 2{Nf - 2)/Nf under the non-anomalous 
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U(l)i^ symmetry. The meson and baryon vevs cannot take arbitrary values but are subject 
to constraints following from (2.1), 

B^^m'} = M^^^Bkt] = B^^^B^^' = B.^Bk,] = m\m'\ - B'^Bu = 0, (2.2) 

where the square brackets denote antisymmetrization. These constraints imply that M, 
and B all have rank less than or equal to 2 and, up to flavor rotations, take the form 




with mim2 = bb and the (Nf — 2) x (Nf — 2) matrix of zeros. For Nf > 3, the quantum 
moduli space is also the same as the classical one [5]. 

For SU(2) superQCD with fundamental flavors M, B, B, and S are thought to generate 
all non-trivial local gauge-invariant operators in the chiral ring of the classical theory [10, 11]. 
When iV/ = 4 or 5 the theory is strongly coupled, and has new massless degrees of freedom at 
the origin of moduli space, so the chiral ring might be deformed or enlarged from the classical 
answer. But for Nf > 6, where the theory is IR free, the classical description is as accurate 
as we like (in the vicinity of the origin of the moduli space) . So we will make the assumption 
that we can write our effective superpotential in terms of just S, M, B, and B. 

However, for Nf > 4, the global symmetries allow infinitely many terms in the effective 
superpotential, making it hard to guess its correct form. So, instead, we use the generalized 
Konishi anomaly (GKA) equations to derive the effective superpotential. If F^{^,Wa) are 
holomorphic functions transforming in the same representation of the gauge group as a chiral 
superfield {i is a flavor index and r an index for the gauge representation), then the GKA 
equation [10] is 

which can be interpreted as the anomalous Ward identity coming from the field transformation 
(5$^ = F.I.. Here Wtrcc is the classical superpotential. The GKA equation is perturbatively one- 
loop exact [10]. It has also been shown [13] that it does not get non-perturbative corrections 
for a U(A^) gauge theory with matter in the adjoint representation as well as for Sp(A^) and 
SO(A^) gauge theories with matter in symmetric or antisymmetric representations. For the 
theories we are discussing here, its non-perturbative status is not known; however, as we 
show below, there is strong evidence, at least for SU(2), that the GKA equations are actually 
non-perturbatively exact. 

Consider now SU(2) superQCD with the classical superpotential 

Wtree = m'f{M{ - M{) + bij{B'^ - +V^Bij - Bij). (2.5) 

Here m^j,bij and ¥^ are Lagrange multipliers constraining the operators M-'-, B^^ and Bij to 
have M-'-, B^^ and Bij as their vacuum expectation values, respectively. (Whenever we need 
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to distinguish an operator from its vev, we put a hat on the operator.) We are looking for 
the effective superpotential WeS as a function of the vevs S, M, B, and B. It follows from 
(2.5) and the nature of the Legendre transform [14, 15, 1] that 



h'^ = - 



lOWeff 



aM^. ' 2 ' - 2 dBij 



(2.6) 



where the factors of 2 come from the antisymmetry of the baryons. 

We now use the GKA equations to determine the dependence of the Lagrange multipliers 
on M, B, B, and S. First set = = Ql in (2.4) yielding 



Mm = S + 2Bb, 



(2.7) 



where we are using a matrix notation on the flavor indices (so that, e.g., the last equation 
Sd^j + 2B'^^bkj). A similar equation, 



stands for M'^j^m^j 



mM = S + 2bB, 



(2.8) 



follows from taking i^^ = = Q^. Two more independent equations follow from taking 
Fr = GabQi and = Qa' and from taking = e"''Q^ and = Q'l, giving 



Bm = -2M^b, 



mB 



-2bM^ 



(2.9) 



We will carry out subsequent calculations at a generic point on the configuration space 
of S, Af, B and B where they are all invertiblc matrices. Note, however, that when Nf is 
odd, B and B, being odd rank antisymmetric matrices, are never invertible. We get around 
this problem by restricting ourselves to an even number of flavors only. Once we have found 
the superpotential for even Nfs, we can add a mass term for one flavor and integrate it out 
to get the effective superpotential for the odd Nf — 1 flavors. 

Therefore, multiplying (2.7-2.9) by appropriate inverses and substituting for the Lagrange 
multipliers m, b, and b using (2.6) gives a set of partial differential equations for WeS 



eflf 



dB^^ 



dB, 



eff 



= S 
= S 

= S 



Integrate the first equation in (2.10) to find 

S 



(BM-^B + M^y^BM-^ 
M-^B{BM-^B + M^)-^ 

M-^B{BM-^B + M^y^BM-^ - M'^ 



Weff = -2 Indet (l + M-^BM-^B^ + G{M, B, S), 



(2.10) 



(2.11) 



where M ^ = {M'^) ^, I is the Nf x Nf identity matrix, and G is an undetermined integration 
function. Comparing the second equation in (2.10) with the derivative of (2.11) with respect 
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to B gives dG/dB = 0. Also, comparing the derivative of (2.11) with respect to M to the 
third equation in (2.10) gives dG/dM = -SM'^, so that 



G = -Slndet{M/A'^) + H{S), (2.12) 

for some undetermined function H{S). The A-dependence was determined by dimensional 
analysis, where A is the strong-coupling scale of the SU(2) superQCD. 

Equivalently, the global flavor symmetry implies that W^s = Wcff (AT, dct M, S) where 
X := M-'^BM-^B. Plugging this functional form into (2.6-2.9) gives simple matrix differ- 
ential equations leading to (2.11-2.12). 

H{S) is determined up to a constant by the U{1)r symmetry. Since R{Wes) = 2, H must 
be linear in S, plus a logarithmic piece to cancel the U(l)ij transformation of the — S'lndet M 
term, giving 

H{S) = (2 - Nf)S[a - ln(5/A3)] (2.13) 

for some undetermined constant a. We can determine a by matching to the Veneziano- 
Yankielowicz superpotential [16], Wyy{S) = 2S[1 — ln{S/AYy^)], for pure SU(2) superYang- 
Mills. It is a short exercise to integrate out the mesons and baryons in (2.11) and match 
strong coupling scales to find a = 1. We therefore find that the effective superpotential is 



Weff = -|ln 



(det M)2 det(I + M-'^BM-^B)] +{2-Nf)S{l-lnS) + {6-Nf)S\nA. (2.14) 



Since S is massive we can integrate it out by solving its equation of motion, dWes/dS 
to find 



Wes{M,B,B) = {2-Nf) 



A^f-^ det My^det(I + M-^BM-^B) 



l/{Nf-2) 



(2.15) 



This superpotential reproduces all known low energy aspects of SU(2) superQCD. The 
easiest way to see this is to convert it to a description which makes the full global flavor 
symmetry manifest. As we mentioned earlier, SU(2) superQCD with Nf fundamental 
and Nf anti- fundamental Qf can be equivalently described in terms of 2Nf doublets Q^, 

N f-\-i ~ h " 

I = l,---,2Nf with Ql = Ql and Qa'" = eabQ-- Hence M, B and B are combined into an 
antisymmetric 2Nf x 2Nf matrix F^-^ = e"''Q^Q^, 

After a bit of algebra^ it is seen that our singular effective superpotential (2.15) can be written 
in terms of this new variable as 

Weff = (2 - Nf) (A^/-6\/ditF) V(Ar/-2) ^ 



^Writo V = (o ^)(y i) with X ■- B-^M, y ■- -B'^M'^. Use the identity det (^ I) = det(l - xy), so 
detV = detBdetSdet(l - xy) = {det Mf det{-y'^x'^ + 1), which gives (2.15). 
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making the SU(2iVj) global symmetry manifest. Indeed, (2.17) coincides with the singular 
effective superpotential found in [7] , and so it satisfies all the checks discussed there: it gives 
rise to the correct moduli space, is consistent under integrating out flavors, and reproduces 
all the higher-derivative F-terms found in [8]. 

The success of this calculation can be taken as evidence that the GKA equations are 
non-perturbatively exact for SU(2) superQCD. 



3. Weff for Nf = Nc + 2 superQCD: non-integrability of GKA equations 

SU(A^c) superQCD has Nf massless quark chiral fields and Nf massless anti-quark chiral 
fields transforming in the fundamental and anti-fundamental representations, respectively. 
Here i = 1, . . . , Nf is the flavor index and a = 1, . . . , Nc is the color index. When Nf = Nc + 2 
the classical moduli space is parameterized by the gauge- invariant vevs of the glueball, meson, 
and baryons defined by 



s 






M] 


•~ QaQji 




Bij 


— -^f -u I. f°-^-°-Ncnki . 


■ ■ WaNc 


B 


Jjki---kNcf n«i 


■ ■ QT' 



The global symmetry of the theory is SU(iV/) x SU(iV/) x 11(1)^ x U(1)r. The 11(1)/? charges 
are R{S) = 2, R{M) = A/Nf, and R{B) = R{B) = 2Nc/Nf. The classical moduh space is 
described by the constraints that M, B, and B satisfy by virtue of their definitions, 

BikM^ = mIb''^ = B[ijBk]i = B^'^B''^^ = B'^Bke - M'^Im-^^} det M = 0. (3.2) 

Square brackets denote antisymmetrization; antisymmetrization on n indices consists of n! 
terms (i.e., with out a factor of 1/n!). 

They imply that by appropriate flavor rotations M, B, and B can be put in the form 




where m is an Nc x Nc matrix and 6, h are numbers satisfying bb = det(m). The classical and 
the quantum moduli spaces are the same [5] , but at the origin there are extra light degrees of 
freedom. At points away from the origin, the only light degrees of freedom are components of 

M, and B. At the origin, SU(A'c) supersymmetric QCD with Nf = A^c + 2 is an interacting 
superconformal field theory for Nc = 2 and 3. For Nc > 4 it becomes strongly coupled, but 
has an IR free dual description in the IR [6]. 
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In this section we will try to construct an effective superpotential in terms of these fields 
which correctly describes the moduli space of vacua for points away from the origin. As in the 
last section, we will use the GKA equations to systematically derive Weff- In fact, the GKA 
equations were used in [9] to construct a set of coupled partial differential equations for the 
effective superpotentials of SU(A''c) supersymmctric QCD. They have been integrated [9] for 
Nf = Nc and Nf = iV^+l where the results are in agreement with those in [5]. Unfortunately, 
the GKA equations are quite complicated for Nf > Nc + 2 flavors. We will show below how 
to simplify the GKA equations when Nf = Nc + 2. 

We briefly recap the derivation of the equations for WeS from the GKA equations [9]. 
The strategy is the same as in the SU(2) case discussed in the last section: start with the 
tree level superpotential 

Wtree = m}(M/ - Mf) + b'^ {Bij - B^j) +\j{Bii - (3.4) 

where 

m^--^ f3 5) 

are Lagrange multipliers enforcing that M- , Bij and Bij be the vevs of the meson, baryon, 
and anti-baryon operators, respectively. There is no need to introduce a Lagrange multiplier 
for S because we are considering points away from the origin and S is massive for these 
points. We get two relations among the Lagrange multipliers by taking Fr to be the quark or 
antiquark field in the GKA equation (2.4), and two more by taking it to be proportional to 
gijfe^2 -^Arc£^^^ ^^^Qa2 . . . and similarly with the Q's. The resulting GKA equations are 

Mlm'] = iS + b''^Bek)S} - 2b"'Bkj, 
miM^ = {S + bu&^)5i - W^Bkj, 



m^^B^^^ = detM, 



m\^Bk^^ = 2hhgM-^lM~^lM-^\ det M. (3.6) 

Note that the right sides of the last two equations, though they are written using M~^, are 
actually polynomial in M. 

Now, the global flavor symmetry implies'^ that 

W.-ii = S fix. det M), where X := — — — , (3.7) 
' " detM ' ^ ' 



^This symmetry argument is not entirely straightforward. The U(1)_b baryon number symmetry imphes 
that for each B there must be an accompanying B in each term. Since Weff is an SU(A//) x SU(Arj) singlet, all 
the flavor indices must be contracted in each term. Contractions with the totally antisymmetric epsilon tensors 
can always be reduced to products of det M and Pf B • Pf B. The only other way to contract indices of B and B 
is with an M as BMB (or its transpose), and since these in turn must be contracted, another factor of M must 
be included. There axe four ways of doing this — BMB and its three cyclic permutations — but upon maicing 
a flavor singlet expression a trace must be taken, so the cyclic order does not matter. Finally, the product of 



Pfaffians of baryons is not independent of X and dot A/, since Pf B ■ Pf B = ^dct(Af^BA/B)/ dot Af. 

Alternatively, one can derive this directly from the GKA equations. Use them to deduce (3.9) and a similar 
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and we are using matrix notation for the meson and baryon fields. The first two equations in 
(3.6) then imply 

>Veff = S W{X) - S In det M + 2S(ln S-1), (3.8) 

where W is to be determined. 

The GKA equations (3.6) imply a matrix differential equation for W{X) as follows. Con- 
tract the ijki indices in the last equation in (3.6) with three M's, giving 2(2 — Nf)b det M = 
tr{mM)M'^BM + M^m^M'^BM-M'^BMmM. Substitute for Mm using the first equation 
to get 

2 6 det M = 



AT — 2 



M'^BM - -^M'^BbBM. (3.9) 



Derivatives of (3.8) with respect to B and B together with (3.5) imply 

2 6 det M = 5* {MBWxM'^ + MW'^BM^), 

2b det M = 5 {WxM^BM + M^BMW^), (3.10) 

where Wx '■= dW/dX. Work at a generic point in parameter space where M, B, B and X 
are all invertible matrices. This is only possible if Nf, and therefore Nc, is taken to be even, 
as in the discussion in the previous section. 

Substitute (3.10) in (3.9) and multiply on both left and right by B to obtain 

BG{X) = -G^{X)B, (3.11) 

where 



G{X) := WxX + ^ 



N —2 
1 + \ tiiWxX) 



X + -^XWxX, (3.12) 



and we have used that BX is antisymmetric. On the other hand, from the definition of X 
(3.7) it follows that BX = X^B, which imphes 

BG{X)=G^{X)B (3.13) 

since G is a function of X alone. (3.13) and (3.11) imply G{X) = 0, which, after being 
multiplied from the right hy X~^, reads 

I / N — 2 \ 2 

Wx + -U + tT{WxX)\ I + j^XWx = 0, (3.14) 

where I is the Nf x Nj identity matrix. This is the matrix differential equation for W{X). 

relation for b, multiply these by B and B, respectively, then substitute the second into the first. One finds 
that bB depends on B and B only through X. Since b ~ dyVes/dB, it follows that the dependence of yVeS on 
B and B is solely through X. 

Note that this symmetry argument is no longer effective when Nf > iV^ + 2. For then B and B have more 
than two indices, and the analog of X is no longer a matrix, but has Nf — Nc — 1 upper and Nf — Nc — 1 
lower (antisymmctrizod) indices. Those objects can bo contractcxl in many incquivalont ways to make flavor 
singlets. This is the source of the difficulty in integrating the GKA equations for Nf > Nc + 2. 
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The trace of (3.14), ti{WxX) = -[2ti(Wx)+Nc+2]/Nc, allows us to eliminate ti{WxX) 
from (3.14), giving N^Wx = {N^ + 2X)-'^[{Nc - 2)tT{Wx) - 2]. The trace of this equation 
allows us to eliminate tr(Wx) in turn, giving the following differential equation for W{X): 

W 2(iV, + 2X)-i 

(7Ve-2)tr((7V, + 2X)-i)-7V,- 

When Nc + 2, we define the matrix Y := (iV^ - 2){Nc + 2X)-\ and substitute it into (3.15) 
to obtain 

dW 1 

where wc have explicitly written the indices to avoid any confusion. This differential equation 
is not intcgrable, as it is easy to check that d'^W/dY/'dYj / d'^W/dYjdY/' . This shows that 
the GKA equations for Weff are not integrable for even values of A^c > 2. 

3.1 Comparison to the SU(2) solution. 

For Nc = 2 we are integrating the same GKA equations as we did in section 2, though in 
terms of the baryons and anti-baryons B and B instead of their Hodge duals B and B. Thus 
the GKA equations must be integrable in this case, and, of course, give the same answer we 
found in section 2, namely, equation (2.15) with Nf = 4. However, there is a subtlety in 
comparing these two computations, which wc will now explain. It will play an important part 
in our discussion of the results of integrating the Seiberg dual GKA equations in the next 
section. 

For Nc = 2, (3.14) is indeed integrable, and integrates to give W{X) = — ^trln(l + X). 
Integrating out S from (3.8) then gives 

A>Veff,su(2) = -2\/detMdet(l + A:)^/^ (3.17) 

where A is the strong coupling scale of the gauge theory. On the other hand, the SU(2) 
effective superpotential for Nf = 4 found in section 2 is 

AyVeff,su{2) = -2VdetMdet(l + M~^BM-'^B)^^'^. (3.18) 

Since the SU(2) baryon fields B and B as defined in (2.1) are Hodge-dual to the S\J (Nc) 
baryons B, B defined in (3.1), and since rank(M) = Nf = A, it follows that 

M-'^BM-^B = ^ tr(A:) - X, (3.19) 
so that the effective superpotential reads 

AWeff,su{2) = -2v'detMdet(l + ^ tr(X) - X)V4. (3.20) 
Apparently the two answers, (3.17) and (3.20), do not agree for general X. 
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The resolution of this paradox is that X is not a general 4x4 complex matrix, but 
satisfies some constraints by virtue of its definition (3.7). In particular, X can be thought of 
as the product of two antisymmetric matrices {M'^BM and B). Such a matrix, though not 
necessarily either symmetric or antisymmetric, has only half the degrees of freedom that a 
general matrix, Y, of the same rank would have. 

To see this, recall that an appropriate similarity transformation G~^YG with G G 
GL(2n, C) will diagonalize Y. Thus the 2n eigenvalues of the general rank 2n matrix Y 
are all its GL(2n, C)-invariants. More properly, a basis of 2n independent symmetric poly- 
nomials of these eigenvalues generates all invariants. This basis can conveniently be taken to 
be tr(yP) for p = 1, . . . , 2n. If, on the other hand, X = AB is the product of two antisym- 
metric matrices A and B, then a GL(2n,C) transformation acts as G~^XG = G^^ABG = 
G~^ AG~^ G^ BG . We can always choose a G = Go so that GqBGq = J where J := I„C5i(j2 is 
the "unit" skew-diagonal matrix. This condition does not fix G, since if H e Sp(2ra, C) (i.e., 
H e GL(2n,C) satisfies H^JH = J), then G = GqH will also give G'^BG = J. For such 
G we have G~^XG = H~^A'H~'^J where A' is the antisymmetric matrix A' = Gq^AGq'^ . 
Now, a Gram-Schmidt othogonalization argument but with respect to the skew product de- 
fined by J shows that H can always be chosen to bring H~^A'H^^ to skew diagonal form 
H~^A'H~^ = diagjai, . . . , a„} ® ia2- Thus X has only n independent (double) eigenvalues, 
and a basis of generators of the GL(2n, C)-invariants of X can be taken to be tr{XP) for 
p = 1, . . . , n. This is half the number of independent invariants of the general rank 2n matrix 
Y, and implies, in particular, that tr(X'') for p > n satisfy additional identities allowing them 
to be expressed in terms of products of traces with p < n. 

For example, for Nf = 4 (ra = 2) the new independent cubic and quartic identities are 
easily found to be 

= StrX^ - etrX^trX + (trX)^ 

= 8trX''-4trX2trX-2(trX2)2-htrX2(trX)2. (3.21) 

When Nf = 4, the identities (3.21) are easily checked to imply that ti{XP) = tr[X'P) for 
p = 1, ... ,4, where X' := ^ tr{X) — X. Thus all invariants made from X' are the same as 
those for X, and in particular, the two forms of the effective superpotential (3.17) and (3.20) 
are equivalent. 

4. Weflf for Nf = Nc 2 superQCD: Seiberg dual analysis 

The IR-equivalent Seiberg dual description [6] of SU(A''c) superQCD with Nf = Nc + 2 is 
SU(2) superQCD with Nf (dual) quarks qf in the fundamental and Nf (dual) anti-quarks 
in the anti-fundamental, and a set of gauge singlets M^j coupled through the superpotential 

W = qlq^Mi, (4.1) 

where i = 1, . . . , Nf and a = 1, 2 are flavor and color indices, respectively. This superpotential 
breaks the global symmetry of the dual theory down to S\J{Nf) x S\J{Nf) x U(l)s x U(l)ij. 
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The dual meson, baryons, and glueball are defined to be 

N]:=qta<lh ^ij-=^amU], B^^ :=e''%qi, S := tv {w^w a)/ {327T^). (4.2) 

This dual theory is IR free when Nf > 6 {Nc > 4), and there are free quarks, anti-quarks and 
gluons at the origin of the moduli space. The U(l)i? charges are R{S) = 2, R{M.) = 4/Nf, 
and R{M) = R{B) = R{B) = 2NJNf. The chiral ring of the dual theory {S,M,B,B) is 
related to that of the direct theory [S, M, B, B) by [1] 

S = -S, M = nM, B = iii-^A^f-^B, B = in-^A^f-^B, (4.3) 
where /x is a matching scale defined by 

A6-JV/ =^^/A6-2iv,^ (4.4) 

with A and A being the dynamical scales of the direct and dual theories, respectively. 
The gauge-invariant form of the classical F-term equations are 

M = BM= MB = 0, (4.5) 

where we are using a matrix notation for the fields. The D-term equations give 

B^^^ = N ^B = B^B = B^B = B®B-^^ ^^^ = (4.6) 

which are the same as the constraints (2.2) of SU(2) superQCD discussed in section 2 with 
the substitution M — A/". The space of solutions to the classical constraints (4.5-4.6) has 
two branches: either M = B = B = and rank(AI) > Nc, or, up to flavor rotations, 

where m is an Nc x Nc matrix and 66 = 0. 

This classical moduli space is clearly not the same as the moduli space (3.2-3.3) of the 
direct SU(A^c) theory. However, the dual theory classical constraints (4.5-4.7) are expected to 
be modified quantum mechanically by strong coupling effects. Because even though the dual 
theory is free at the origin of moduli space, arbitrarily small vcvs for Ai destabilize the free 
fixed point by giving masses to the dual quarks through the superpotential coupling (4.1), so 
the theory flows to strong coupling for large enough rank of M . Strong coupling effects are 
argued in [6] to generate the constraints (3.2) of the direct theory. 

In particular, the branch of (4.5-4.7) with rank(AI) > Nc is lifted, and the 66 = 
constraint on the other branch is deformed to 66 = det(m). As an example — which will be 
useful later — of how the classical constraints are modified by strong coupling effects, note 
that if the singlet vev M. is given a generic value {e.g., by constraining it with a Lagrange 
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multiplier term as we will do below) with B = B = 0, the superpotential (4.1) gives mass 
to all the dual quarks. The dual theory thus flows to SU(2) super Yang-Mills in the IR with 
glueball vev S = (A^~^f detM)^^^ generated by a Veneziano-Yankielowicz superpotential 



eff,VY = 



where A is the strong-coupling scale of the dual theory. The scale of the glueball appearing 
this superpotential is determined by one loop matching Ay^ = A^^^-' detA^. Integrat- 
ing S out of (4.8) then generates an effective superpotential for the singlet field given by 
2A(6-^/)/2VdetA^, thus lifting the rank(A^) = Nj region of the classical moduli space. 

4.1 Derivation of WeS 

We now wish to find an exact effective superpotential for the dual theory which reproduces 
the strong quantum effects described in the last two paragraphs. As in previous sections we 
start with a tree level superpotential 

Wtree = M^M^ + w^^iMi - M^) + b'^ {Bij - Bij) + hij{B^i - B'^), (4.9) 

where 



, _ dWeS ,ij iSWeff r l9Weff ..... 



are the by-now familiar Lagrange multipliers. Wc have only specified the vevs of i3, and 
i3, because they completely parametrize the moduli space. 

Using the GKA equations (2.4) in precisely the same way as before, we find a set of 
equations similar to those, (2.7-2.9), found in section 2, 

MN = 5 + 26i3, NM = 5 + 2i36, BM^ = -2J\fb, M^B = -2bJ\f, (4.11) 

together with the equation coming from the variation of the singlet field M: 

M = -m. (4.12) 

Eliminating J\f by plugging (4.12) into (4.11) gives a set of partial differential equations for 
the effective superpotential 

- Mm = S + 2bB, -mM = S + 2Bb, BM^ = 2mb, M^B = 2bm. (4.13) 

To solve for the effective superpotential we multiply the third equation in (4.13) from 
the right by B and from the left by M to obtain MBM^B = 2MmbB = -2(51 + 2bB)bB, 
where in the second equality we used the first equation in (4.13) to eliminate Mm. Solving 
this quadratic equation for bB, we find 



bB 
Bb 



= ^ (^-S + - iMBM^B^ , Mm = ^ (^-S - Vs^ - 4MBM'^B^ , 

= ^ (-5 + - 4BM'^BM^ , mM = ^ (-S - Vs^ - 4BM^BM^ , (4.14) 
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where the second Une comes from solving similar quadratic equations for Bb and rriAi. 

The matrix square roots in the above expressions need some explanation. In order to 
make sense of them, consider a region in the configuration space where the magnitude of 
each element in the matrix MBAd^'^B is much smaller than S^. We can then expand the 
the square root as a power series, V S'^ - AMBM^B = S\/l{l - 2S^'^MBM^B + •••), and 
by analytic continuation we extend the result to include all points on the parameter space. 
However, for this to be a definition of the square root, we still need to determine VI. In 
general \/i = 2P — I where P can be any projection matrix (P^ = P). Which P should 
we use? The following argument shows that we have to take \/l = ±1. As we saw in the 
discussion surrounding (4.8), we expect to generate a non-zero S at points in the parameter 
space where ^A ^ Q but B = B = so this is a suitable region to evaluate the square roots. 
At such points (4.13) implies that mM = Mm = —S and BM^ = M^B = 0. But this is 
only consistent with (4.14) if ^/SH. = SI implying that \/l = ±1. 

It is straightforward to integrate (4.14) for the effective superpotential to get 



eff 



In det 



1 



AMBM^B - S \ _^^^^^^^^-ej^jsj^Tjg^^^^^^^^^^-ij^^ 

AMBM^B + S 4 



1 



+ - tr - AMBM^B + ^(4 - Nf)S a - ln(5/A^ 



(4.15) 



We used the 11(1)^ symmetry to fix the S dependence up to an undetermined integration 
constant a. Evaluating WeS at B = B = gives, after a somewhat delicate cancellation. 



>Veff(H = B = 0) = 25 a + (1 



a)—r — In 
4 



(4.16) 



Comparing to the answer (4.8) expected from the strong coupling analysis, fixes a = 1. 
Note that this limiting form (4.16) is already a check that the effective superpotential (4.15) 
is consistent with the quantum modified constraints. In particular the S In dot M term in 
(4.16) serves to lift the whole rankM = Nf region of the classical moduli space (4.5-4.7), in 
accordance with the expected quantum constraints (3.2-3.3). 

So, our final result for the effective superpotential for the Seiberg dual theory is 



Weflf = ^ det 



AMBM'^B - S 



VS^ -AMBM'^B + S^ 



- - In det (mBM'^B^ + 5 In det M 



1 C 

+ - tr - AMBM'^B + (4 - iV/)-(l - In 5) + (6 - Nf)S In A. (4.17) 

Using the chiral ring mappings (4.3), this can be expressed in terms of the fields of the direct 
theory. Note that since S = —S, the definition of the branch of the square root made above, 
VS^I = SI, now becomes V S"^! = —SI. We make this minus sign explicit by changing the 
signs of all square roots and keeping the convention VSH = +SI. We then find 



In 



det 



\/ s'^ + Ax + s\ det X 



Vs2 + Ax -s det^ M 



^ tr +Ax + {Nf - 4) I [1 - In s] , (4.18) 
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where we have defined the shorthands 



A^f-^WeS, s:=A^f-^S, x:=MBM^B, (4.19) 



to avoid having to write factors of A. This is the effective superpotential for SU(A^c) superQCD 
with Nf = Nc + 2. 

4.2 Integrating out the glueball field 

Away from the origin of moduh space, the glueball s is expected to be massive. Solving its 
equation of motion, dwes/ds = 0, gives s = s^{M,B,B) where is defined implicitly by 



sr^-'^ = det ( ^^S±il ) (4.20) 

det^M ^ ' 

Substituting this in gives the effective superpotential as a function of meson and baryons 
only, 

(4.21) 



1 + ^tr (Jl + A^-I 



4 \y si 



The relation (4.20) gives s^, as a complicated function of M, B, and B, which makes it difficult 
to deduce the equations of motion for these fields from WcS- 

Before solving (4.20) for s*, we can extract some of the constraints that M, B, and B 
must satisfy on the moduli space. Since WeS is singular, as shown in [7] it must be regularized 
in order to extract its physical predictions. The idea behind the regularization is to deform 
Weff by introducing some regularizing parameters /x, ^ as follows 

WeS ^ = + I^Mj + p'^Bij + Pij&K (4.22) 

Generic small values of ^, /3, and /?, will fix M, B, and B to some values which will be shifted 
from the moduli space (the extrema of tfcfr) by positive powers of the regularizing parameters. 
Thus as /3, /3 ^ 0, {M, B, B} will approach some point on the moduli space. However, the 
specific point reached on the moduli space will depend on how the regularizing parameters 
scale to zero: as different orders of limits of /x, /?, /3 — >^ are taken, the whole moduli space 
will be scanned. (Note that some vanishing limits of the rcgularizcrs can also send M, S, or 
B to infinity; this is unavoidable since the moduli space itself stretches off to infinity.) 

This regularization procedure should be compared with the trick wc used of introducing 
Lagrange multipliers (4.9) to derive equations for w^^ from the GKA equations. With the re- 
placement of the Lagrange multipliers with the regularizing parameters, {m, 6, 6} {//, /3, /?}, 
the identification of the Lagrange multipliers as derivatives of WeS-, (4.10), is now interpreted 
as the regularized equations of motion. Thus, extremizing w^^'^ using (4.18) gives 



4l3B = s - -x/s^ + 4x, 2Mfi = s + Vs^+ix, 

4Bp = s- ^s'^ + 4j/, 2ijlM = s + ^J + Ay, (4.23) 
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where y := M^^xM. (4.23) is equivalent to (4.14), the original differential equations — with 
the operator mappings (4.3) — we integrated to get in the first place. Then, reversing the 
manipulations which led from (4.13) to (4.14) gives 

Mn = s-2/3B, nM = s-2Bp, BM^ = -2fiP, M^B = (4.24) 

These show that independent of how IJ,,P,P —>■ 0, we always end up with BM'^ = M'^B = 
0. This implies in particular that x = —4Pfi^P^, and so also always vanishes with the 
regularizing parameters. Also, if as //,/?, /3 — M, B, and B remain finite, then s — as 
well. 

This shows that the extrema of the superpotcntial (4.18) satisfy the constraints BM'^ = 
M^B = s = 0, which are, indeed, part of the constraints (3.2) describing the Nf = Nc + 2 
superQCD moduli space. The remaining constraints should also follow from the effective 
superpotential. However, they are much harder to derive, as they require solving for s* in 
(4.20). Now the argument of the last paragraph shows that s* = 0, but, because of the 
singular nature of the superpotential, it is incorrect to simply plug this value into (4.21) to 
find tt^cff. Instead, should be found for generic M, B, and B, and then the extrema of the 
resulting effective superpotential can be analyzed by regularizing it, as above. 

The equation (4.20) can be solved systematically for by assuming that (all the eigenval- 
ues of) ( := xjsl <^ 1 and expanding in powers of this parameter. The leading order solution 
is si = detM, which can be checked to be consistent with the assumption that ||^|| ^ 1. 
Change variables to 

so that (4.20) and (4.21) become 

= det X det ( :^^^±A = 3-2^/ det (V^^T^X + a? , (4.26) 
\ A/fT^ -I- AX — rr \ / 



= -2Vdet M 



C7 + ^ tr (^V<J^ + 4X - a) 



(4.27) 



Expand the right side of (4.26) in powers of X/a and solve it consistently order-by-order in 
a power series expansion fi^ = 1 + • • • to get 

2 ^ ^ tvX (trX)2 3tr(X2) (trX)^ 3trXtr(X2) 5tr(X3) 

2 8 4 12 4 3 ^ ' ' 

_ 9(trX)^ 27(trX)2tr(X2) _ 27ti{X^f _ 5trXtr(X3) 35tr(X^) ,^.^5x 
128 32 32 2 8 ^ ^' 

Plugging this into (4.27) gives 



Wes = -2Vdet M 



tiX (trX)2 _ tr(X2) 5(trX)3 _ 3trXtr(X2) tr{X^) 
^ ~T~ ~^ 32 8 ^ 384 32 ^ 6 

^ 49(trX)^ _ 21(trX)2tr(X2) ^ 9tr(X2)2 ^ 5trXtr(X3) 5tr{X^) 



6144 256 128 24 16 

+ 0{X^)\. (4.29) 
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This is the answer, to order X'^, for the effective superpotential for Nf = Nc + 2 superQCD 
found by integrating the Seiberg dual GKA equations and integrating out the glueball. 

Alternatively, we can derive a differential equation satisfied by w^s as a result of in- 
tegrating out s. Two of the equations of motion that we originally integrated to get the 
superpotential (4.23) become, when written in terms of a and X, 

Af3B = VdetM (a - ^ a"^ + 4x) , 

2M/X = VdetM (g + \/ a"^ + 4x) . (4.30) 

On the other hand, upon integrating out s, we have seen that w^^ takes the form 

u;eff = VdetM /(X) (4.31) 

for some function /. (This form just follows from the symmetries.) Two of the equations of 
motion following from this form of lUeff are 

4/?S = VdetM ( 4X/' ) , 

2M\i = VdetM (-4X/' - / + 2 tr(X/')) , (4.32) 

where /' is the matrix derivative df /dX. 

Equating and adding (4.30) and (4.32) gives 

a = -i/ + tr(X/'), (4.33) 

while equating and multiplying (4.30) and (4.32) gives 

1 = /'(/- 2 tr(X/') + 4X/'). (4.34) 

This is a nonlinear first order (matrix) differential equation for WeS with the glueball integrated 
out. We do not know how to integrate this equation in closed form when rank(X) > 1. But 
it is straightforward to check that a series expansion of the solution to (4.34) with boundary 
condition /(O) = —2 reproduces (4.29). 

It may be clarifying to note that (4.34) has a one-parameter family of solutions. For if 
f{X) is a solution, then so is 

/(„)(X):=a/(a-2x) (4.35) 

for any a G C*. By (4.33) this implies that a{X) changes to aa{a~'^X); but it is easy to check 
that the a equation of motion (4.26) does not have this symmetry. The algebraic equation 
(4.26) determining a has more information than the differential equation (4.34), and so picks 
out a single instance of the family of solutions (4.35), namely the one obeying the boundary 
condition /(O) = —2. (For example, — 2tr(V2X) solves (4.34) but does not satisfy (4.26), so 
is not a physical solution.) 
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4.3 Comparing to the direct result when Nf = 4 

Is (4.29) the correct superpotential? A basic check is to see whether its extrema (computed 
by appropriately regularizing, as explained above) reproduce the moduli space given by the 
constraints (3.2). This seems a very difficult check to perform since we do not have a closed 
analytic form for WeS- 

However, there is one case where we can carry out a non-trivial check. When Nf = 4, 
Nc = 2, so in this case we should reproduce the superpotential (3.17) found in sections 2 and 
3 for the SU(2) theory. Expanding (3.17) in powers of X, wc find 

AWe«,su(.)=-2VditM[l + ^ + 4^-^ + i^ 32^ + ^ 

(trX)4 _ (trX)2tr(X2) ti{X^)^ tiXti{X^) _ tiiX^) 
6144 256 ^ 128 48 16~ 

+ 0{X^)\. (4.36) 

Though it does not coincide with the expansion (4.29) starting at order X^, we must bear 
in mind the identities (3.21) that traces of powers of X satisfy starting at cubic order. One 
finds that the difference between (4.36) and (4.29) is proportional to these identities, and 
so vanishes. Thus the effective superpotential found by integrating the dual GKA equations 
matches the correct result at least to quartic order in X. 
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